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Abstract 

We calculate accurate bound states and resonances of two interest- 
ing perturbed Coulomb models by means of the Riccati-Pade method. 
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H , This approach is based on a rational approximation to a modified 

logarithmic derivative of the eigenfunction and produces sequences of 
roots of Hankel determinants that converge towards the eigenvalues 
of the equation. 



1 Introduction 

In a most interesting series of papers Killingeck et al [1-3] and Killingbeck 
[4, 5] have shown that perturbation theory and the Hill-series method are 
suitable tools for the calculation of bound states and resonances of simple 
quantum-mechanical models. In order to obtain the complex eigenvalues 
that correspond to unstable states they resort to a complex parametrization 
of the methods that they call "complexification " . 

Another approach that proves useful for the accurate calculation of bound 
states and resonances is the Riccati-Pade method (RPM) based on a rational 
approximation to a modified (or regularized) logarithmic derivative of the 
eigenfunction [6-15]. In this paper we apply the RPM to the interesting 
perturbed Coulomb problems discussed recently by Killingeck [5] with the 
purpose of challenging the recently developed asymptotic iteration method 
[16-22]. 

In Section[2]we outline the main features of the RPM. In Section[3]we dis- 
cuss a perturbed Coulomb model with interesting bound states. In Section H] 
we calculate the resonances for a slightly modified model with continuum 
states. Finally, in Section \5\ we draw conclusions on the performance of the 
RPM. 

2 The Riccati-Pade method (RPM) 

Suppose that we want to obtain sufficiently accurate solutions to the eigen- 
value equation 

ij"{x) + Q{x)ij{x) = 0, Q(x) =E- V{x) (1) 



where Q{x) can be expanded as 

oo 

Q{x) = Y: Qj-2X^'-' (2) 

i=o 

about a; = 0. We transform the hnear differential equation ([T]) into a Riccati 
one for the modified logarithmic derivative of the eigenf unction: 

-- ' ^'(^) (3) 



^''^' = x- 


i^ix) 


On substituting ([3]) into ([1]) we obtain 




fix) + -fix) - fix)' - 


-Oix) '^'-^^ 


^\^) ^2 



(4) 

We choose s(s — 1) = —Q^2 in order to remove the pole at origin, and, as a 
result, the expansion 

oo 

fix) = x^-' y: fj^^' (5) 

j=0 

for the solution to the Riccati equation (jl]) converges in a neighbourhood 
of X = 0. Notice that if we substitute the expansions ([2]) and i^ into the 
Riccati equation (jlj) we easily obtain the series coefficients fj as a function 
of E and the known potential parameters Qj. 

We rewrite the partial sums of the expansion ([5]) as rational approxima- 
tions x^~^[N + d/N]iz), where z = x^, in such a way that 

[N + d/N]iz) = ^;\' ^ = E f,z^+Oiz'^^'^') (6) 

In order to satisfy this condition the Hankel determinant Hf,, with matrix 
elements fi+j+d-i, hj = 1, 2, . . . , D, vanishes. Here, D = A^ + 1 = 2, 3, . . . 
is the determinant dimension and d = 0, 1, . . .. The main assumption of the 
Riccati-Pade method (RPM) is that there is a sequence of roots E^^''^ of 



Hf){E) = for D = 2, 3, . . . that converges towards a given eigenvalue of 
equation ([1]). Comparison of sequences with different d values is useful to 
estimate the accuracy of the converged results. 

Notice that we do not have to take the boundary conditions explicitly into 
account in order to apply the RPM; the approach selects them automatically. 
In addition to the answers expected from physical considerations, the RPM 
also yields unwanted solutions as shown below. 

3 Model with bound states 

From the ansatz (f{r, X) = r exp (— r — Ar^) and the equation ip"{r, A)/ [2(f{r, A)] 
V{r, A) — E{\) we derive a potential-energy function V{r, A) = — l/r + 2Ar + 
2A2r2 if E{\) = -1/2 + 3A. For A > y?(r, A) and E{\) are a pair of 
eigenfunction and eigenvalue of the Schrodinger equation with the potential 
V{r, A). For A < E{X) = —1/2 + 3A is not an eigenvalue of the Schrodinger 
equation because the corresponding eigenfunction ip{r, A) is not square in- 
tegrable. Curiously enough, e(A) = —1/2 — 3A is close to the ground-state 
eigenvalue of the Schrodinger equation 

ij"{r) +2[E- Vi{r)] V^(r) = 0, Vi{r) = -- - 2Ar + 2AV^ A > (7) 

when A is sufficiently small. Killingbeck [5] calculated the energy shift A (A) = 
E{X) — e(A) very accurately for several values of A by means of the Hill-series 
method. 

Our interest in this model stems from the fact that 1/r — (p'{r, —A)/ 
ip{r, —A) = 1 — 2Ar is an exact rational function and therefore e(A) will always 
be a root of the Hankel determinants even though it does not correspond to 



a square-integrable eigenfunction if A > 0. This unwanted solution will 
appear as an exact multiple root of the Hankel determinant, very close to 
the physical one when A is close to zero. 

If A < 2/27 the potential-energy function ([7]) exhibits three stationary 
points: a minimum at ri < 0, a maximum at r2 > 4 and a shallow minimum 
at ra > 4. On the other hand, there is only a minimum at ri < when 
A > 2/27. Obviously, only the stationary points at r > make sense from a 
physical point of view, and we expect the RPM to yield better results in the 
latter case. The expansion of the solution to the Riccati equation about r = 
will require many terms in order to take into account the shallow minimum 
that will move away from origin as A decreases. In this case we expect to face 
the necessity of Hankel determinants of greater dimension in order to obtain 
the shift to a given accuracy as A decreases. This unfavourable situation is 
an interesting test for the RPM that has not been applied to this kind of 
problems before. 

The Hankel determinants are polynomial functions of A and A. For ex- 
ample, H^{A, A) = A^^^Pd{A, A) and, therefore, the approximation to the 
energy shift is given by a root of Pd(A, A) = that approaches the multiple 
root A = as A decreases. Table [T] shows A (A) for some values of A calcu- 
lated with determinants of dimension D < 20. In order to estimate the last 
stable digit we compared the sequences of roots with d = and d = 1. As 
expected from the argument above, the accuracy decreases as A decreases if 
we do not increase the maximum value of D consistently, but in all cases we 
have verified that there is a sequence of roots of the Hankel determinants 
that converge towards the ground-state eigenvalue. Present results agree 



with those calculated by Killingbeck by means of the Hill-series method [5]. 

4 Model with no bound states 

It has already been shown that the RPM is a most efficient tool for the 
calculation of resonances in the continuum of simple quantum-mechanical 
models [11,14,15]. However, for completeness in what follows we consider 
the potential-energy function 

V2(r) = — + 2Ar-2AV2 (8) 

r 

that is closely related to the preceding one but does not support bound 
states because it is unbounded from below as r ^ oo. In this case we 
expect unstable or resonant states with complex eigenvalues that correspond 
to tunnelling from the Coulomb well. 

Table [2] shows present results obtained from Hankel sequences with D < 
20. As in the preceding example we compared sequences with d = and 
d = 1 in order to estimate the last stable digit. Our results agree with those 
reported by Killingbeck [5], except for A = 0.08. While the ffist digits of the 
imaginary part of our eigenvalue agree with those in Killingbeck's Table 3 [5], 
the real part is completely different. The disagreement is due to a misprint 
in Killingbeck's Table 3 for that particular entry. In fact, we have found 
that the real part of the eigenvalue reported by Killingbeck for A = 0.08 
corresponds to A = 0.05 instead, as shown in present Table O 



5 Conclusions 

We have shown that the RPM is suitable for the accurate calculation of bound 
states and resonances of perturbed Coulomb problems. The first model, 
equation ([7]), considered in this paper exhibits interesting features that were 
not faced in previous applications of the RMP [6-15]. A shallow minimum 
that moves forward from origin as the potential parameter A decreases makes 
it necessary to resort to Hankel determinants of increasing dimension in order 
to obtain eigenvalues of a given accuracy. On the other hand, we had already 
proved that the RPM is suitable for the calculation of resonances in the 
continuum, and we simply confirmed this strength of the approach by means 
of the second model chosen above. 

The applicability of the RPM is not restricted to eigenvalue equations. We 
have recently applied it to several ordinary nonlinear differential equations 
[23] . Since most of them are not Riccati equations we called this variant of 
the method Pade-Hankel, but the strategy is basically the same outlined in 
this paper. 
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Table 1: Energy shift A (A) for the ground-state energy of the perturbed 
Coulomb model ([7]) 



A A(A) 



0.10 3.41730960373299 10^^ 

0.09 2.31341988422733 10'^ 

0.08 1.4212168993068 10-2 

0.07 7.546639486534 lO^^ 

0.06 3.1738752354 10"^ 

0.05 8.93101948 10"^ 

0.04 1.1819718 10"^ 



Table 2: Resonance for the Is state of the perturbed Coulomb model ([8]) 

A ReE ImE 

0.10 -0.27519233330828482428 1.3918964850900 10"^ 

0.09 -0.29265795893536614770 7.9213310722 10-^° 

0.08 -0.31105186469292522577 2.094858859 10^^^ 

0.05 -0.372260539194895485 



10 



